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We present a measurement of the time-dependent CP- violating asymmetry in K*^j decays 

with K*° K'^tt" based on 232 miUion T{4S) BB decays collected with the BABAR detector 
at the PEP-II asymmetric-energy e'^e~ collider at SLAC. In a sample containing 157 ± 16 signal 
decays, we measure Sk-o^ = -0.21 ± 0.40 ± 0.05 and Ck-o^ = -0.40 ± 0.23 ± 0.03, where the 
first error is statistical and the second systematic. We also explore B" KsT^'^'f decays with 
1.1 < mjfo^o < 1.8 GeV/c^ and find 59 ± 13 signal events with Sko^o-, = 0.9 ± 1.0 ± 0.2 and 
Cjfo^o^ = -1.0±0.5±0.2. 

PACS numbers; 13.25.Hw, 13.25.-k, 14.40.Nd 



The decay transition 6 ^ S7 is sensitive to contribu- 
tions from physics beyond the Standard Model (SM) 0|. 
There has been extensive experimental and theoretical 
investigation of the inclusive decay rate B{B Xsj), 
which to date shows no significant deviation from the 
SM 0. Various new physics scenarios can accommo- 
date large deviations from the SM in other b ^ sj decay 
properties as well, in particular in CP-violating (CPV) 
asymmetries and the polarization of the final state pho- 
ton 0. The photon polarization in 5 ^ 57 (6 ^ 57) 
is predominantly left handed (right handed) in the SM. 
As a consequence, in the exclusive decay B'^ {K^Tr'^)'j 
interference of the amplitude for the direct decay and 
the amplitude for the decay via — B'^ mixing is sup- 
pressed. Therefore, time-dependent CP-violating asym- 
metry is expected to be small 0, Sk-'O-j ~ sin2/3 ~ 
—0.04, where (wf,) is the mass of the s (6) quark, 
(3 = arg(— T4dl/J/VtdVjj) and V is the quark mixing 
matrix Any significant deviation that goes beyond 
possible hadronization corrections of order 0.1 j5|] would 
indicate phenomena beyond the SM. 

In this Communication we report new measurements of 
the time-dependent CPV asymmetry in B'^ —> K 
based on 232 million T(4S') —f BB decays collected with 
the BABAR detector at the PEP-II asymmetric-energy 
e+e" collider at SLAC. Measurements of the CPV asym- 
metry in B^ K*^j, the subset of events with 0.8 < 
it^k^tt" < 1-0, have previously been reported by BABAR 
on 110 fb"^ 1?^ and BELLE on 253 fb"^ 0. The BELLE 
collaboration has also reported a measurement of inclu- 
sive B° i^^7r"7 with 0.6 < m^o^o < 1.8 GeV/c^ Q. 
The latter measurement is motivated by a recent theo- 
retical result that indicates that all contributions to the 
K^TT^j final state have the same CP eigenvalue Q, so 
that beyond-the-SM effects can be discovered even if the 
m^o^o resonance structure is not resolved. Since the 
correctness of such an averaging procedure is still un- 
der discussion we present our results for events with 



an invariant mass of the i^gTr" pair near and above the 
K*{892)'^ resonance separately. For simplicity we refer 
these two contributions as "K*" and "non-i^T*' respec- 
tively. 

The BABAR detector is fully described in Ref. [13. The 
components that are most important for this analysis 
are a five-layer double-sided silicon micro-strip detector 
(SVT), a 40-layer drift chamber (DCH) and a CsI(Tl) 
electromagnetic calorimeter (EMC). For event simulation 
we use the Monte Carlo event generator EVTGEN ^ 
and GEANT4 pj]. 

At the T{AS) resonance time-dependent CPV asymme- 
tries are extracted from the distribution of the difference 
of the proper decay times At = tcp — ^tag, where tcp 
refers to the decay time of the signal B (Bcp) and itag 
to that of the other B (Btag). The At distribution for 
Bcp f follows 



-\At\/TB 

V±{At) = — X 

4tj 

[1 ± 5. 



4tb 

1/ sin {ArridAt) =F C/ cos (Am^At) ] , 



(1) 



*AIso with Universita di Perugia, Dipartimento di Fisica, Perugia, 
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where the upper (lower) sign corresponds to Btag decay- 
ing as B° (B°), tb is the B° lifetime and Am^ is the 
mixing frequency. The cocfhcients C/ and Sf can be ex- 
pressed in terms of the B'^-B'^ mixing amplitude and the 
decay amplitudes for B^ f and B^ ^ f j^^l- Direct 
CP violation in the decay B° — > / results in a non- 
zero value of C/. For B° K*^^ direct CP violation 
is constrained by measurements of the partial rate asym- 
metry in decays with K*^ K^it~ , A^,o^ = —Ck'O-j = 
—0.010 lb 0.028 Il4l . which is in good agreement with the 
SM prediction [l5l| . 

We search for B° K^n^^ decays in BB candi- 
date events, which are selected based on charged par- 
ticle multiplicity and event topology ■ Candidates for 
— > 7r+7r~ are formed from pairs of oppositely charged 
tracks with a vertex probability larger than 0.001, a 
TT+TT" invariant mass 487 < m^+^- < 507 MeV/c^ and 
a reconstructed decay length greater than five times its 
uncertainty. Photon candidates are reconstructed from 
clusters in the EMC that are isolated from any charged 
tracks and have the expected lateral shower shape. We 
form TT^ — ^ 77 candidates with an invariant mass 115 < 
TO-y-y < 155 MeV/c^ and energy B^o > 590 MeV from 
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pairs of candidate photons each of which carries a min- 
imum energy of 30 MeV. For the photon from the B 
decay, the so-called primary photon, we require an en- 
ergy in the e+e~ frame of 1.5 < E* < 3.5 GeV. We 
veto primary photons that form a tt*' ^ 77 (77 — > 77) 
candidate with invariant mass 115 < m^^ < 155 MeV/c^ 
(470 < rrijj < 620 MeV/c^) when combined with another 
photon with energy > 50 MeV {E^ > 250 MeV). 

To identify 5° decays in K'^n'^'y combi- 
nations we use the energy-substituted mass 
TOES = ^A^/^^-^^Pi^~PB)yEf^^]^ and the energy 
difference AE = Eg- \fsl2. Here [Ei, pi) and {Eb,Pb) 
are the four-vectors of the initial e~^e~ system and the B 
candidate, respectively, is the center-of-mass energy, 
and the asterisk denotes the e~^e~ rest frame. For signal 
decays, the toes distribution peaks near the B mass 
with a resolution of about 3.5 MeV/c^ and AE peaks 
near OMeV with a resolution of about 50 MeV. Both 
Toes and AE exhibit a low-side tail from energy leakage 
in the EMC. We require 5.2 < toes < 5.3 GeV/ and 
|A_B| < 250 MeV, which includes the signal region as well 
as a large "sideband" region for background estimation. 
We also require | cos6'|j| < 0.9, where 9*^ is the angle of 
the B candidate with respect to the e~ momentum in 
the e+e~ rest frame. Finally, for the subset of events 
with mp^a^o < 1.1 GeV/c^, we require \ cos6k*\ < 0.9, 
where 6^' is the angle between the K'^ and the primary 
photon in the K^^tt^ rest frame (the "helicity" angle). 

Event topology is exploited to further suppress the 
background from continuum e^e^ ~^ QQ il — d, s, c) 
events. We calculate the ratio L2/LQ of two moments 
defined as Lj = |p*|| cos0*p, where p* is the mo- 
mentum of particle i in the e^e~ rest frame, 9* is the 
angle between p* and the thrust axis of the B candidate 
and the sum runs over all reconstructed particles except 
for the B candidate daughters. We require L2/L0 < 0.55, 
which suppresses the background by more than a factor 3 
at the cost of approximately 10 % signal efficiency. After 
all selections are applied the average candidate multiplic- 
ity in events with at least one candidate is approximately 
1.1. We select the candidate with a reconstructed 7r° 
mass closest to the expected value and if ambiguity per- 
sists we select the candidate with mass closest to the 
expected value. 

Selected events are divided in events with 0.8 < 
TOj^o^o < 1.0 GeV/c^, where signal decays are predom- 
inantly B^ i^*°7, and events with 1.1 < mj^o^a < 
1.8 GeV/c^, where the contribution from K*{892) is 
small. In the data we find respectively 1469 and 2629 
candidate events in these categories. The selection effi- 
ciency for B^ — > K*'^'y, evaluated with simulated events, 
is approximately 16%. Using the current world average 
for the branching fraction [Tg we expect 176 ± 18 signal 
events. Compared to our previous measurement the 
current event selection is more effective in suppressing 
background from B decays, leading to a reduced sys- 
tematic uncertainty from an eventual CPV asymmetry 



in the background without a significant loss in statisti- 
cal sensitivity. The selection efficiency for i?° Kgir'^'^ 
events with 1.1 < mj^a^o < 1.8 GeV/c^ is approximately 
15%, but depends on the helicity structure. Besides the 
K*{892) the only observed Kit resonance in _B ^ Ktt^ 
decays is the if2(1430). Using the world average for 
the B° X2(1430)°7 branching fraction 22] we expect 
24 ± 7 events. However, since upper bounds on other res- 
onances are weak, the actual observed signal yield may 
be appreciably higher. 

For each B candidate we examine the remaining tracks 
in the event to determine the decay vertex position and 
the flavor of Stag ■ Using a neural network based on kine- 
matic and particle identification information jl8l | each 
event is assigned to one of seven mutually exclusive tag- 
ging categories, designed to combine flavor tags with sim- 
ilar performance and At resolution. We parameterize the 
performance of this algorithm in a data sample (i?flav) of 
fully reconstructed B° D''*^~tt+ / p+ /a^ decays. The 
average effective tagging efficiency obtained from this 
sample is Q = Ec4(l " 2u;'=)2 = 0.305 ± 0.004, where 
e| and w"^ are the efficiencies and mistag probabilities, 
respectively, for events tagged in category c = 1, ... 7. 

The proper-time difference is extracted from the sep- 
aration of the Bcp and Stag decay vertices in a manner 
analogous to Ref. The i3tag vertex is reconstructed 
from the remaining charged particles in the event [l^ . To 
reconstruct the Bcp vertex from the single K'^ trajectory 
we exploit the knowledge of the average interaction point 
(IP), which is determined from the spatial distribution of 
vertices in two-track events and is calculated separately 
for each 10-minute period of data-taking. We compute 
At and its uncertainty from a geometric fit to the 
T(45) B^B^ system that takes this IP constraint into 
account. We further improve the At resolution by con- 
straining the sum of the two B decay times (^cp +^tag) to 
be equal to 2 r^o with an uncertainty \/2 t^o . We have 
verified in a Monte-Carlo simulation that this procedure 
provides an unbiased estimate of At. 

The per-event estimate of the uncertainty on At re- 
flects the strong dependence of the At resolution on the 
flight direction and on the number of SVT layers tra- 
versed by the decay daughters. In about 70 % of the 
events both pion tracks are reconstructed from at least 
4 SVT hits, leading to sufficient resolution for the time- 
dependent measurement. The average At resolution in 
these events is about 1.1 ps. For events that fail this cri- 
terion or for which a{At) > 2.5 ps or |At| > 20 ps, the 
At information is not used. However, these events still 
contribute to the measurement of Cn'O^, which can also 
be extracted from flavor-tagging information alone. 

Signal yields and CPV asymmetries are extracted us- 
ing an unbinned maximum- likelihood fit to toes, AE, 
L2/L0, flavor-tag. At and a{At), as in Ref. For 
the analysis of the B° K*'^^ sample m^o^o is also 
used in the fit. Because we expect a contribution from 
other B decays {"BB background"), we allow the fit to 
extract the fraction of such decays as well. We have 



6 



verified using fits to simulated samples that the corre- 
lation between the observables is sufhciently small that 
the event likelihoods for signal Vs, BB background 
and continuum background Vqq can be described by the 
product of onc-dimensional probability density functions 
(PDF). The PDFs for signal events and BB background 
events are parameterized using either the -Bflav sample 
(for the flavor-tag efflciency, mistag probabilities and Ai- 
resolution function) or simulated events. For the contin- 
uum background, we select the functional form of the 
PDFs in background-enhanced samples. We exploit the 
large fraction of background events in the final sample 
to extract the background parameters along with the 
physics measurements in the fit. The asymmetry in the 
rate of i?" versus tags in background events is also 
extracted from in the fit. 

The PDF for the At of signal events and BB back- 
ground events is obtained from the convolution of Eq. ^ 
with a resolution function TZ{5t = Ai — Aitruci cta*)- The 
asymmetries S^^ and Cg^ for the BB background are 
fixed to zero in the fit, but we account for a possible 
deviation from zero in the systematic uncertainty. The 
resolution function is parameterized as the sum of three 
Gaussian distributions jl6| |. The first two Gaussian dis- 
tributions have a width proportional to the reconstructed 
a At and a non-zero mean proportional to a At to account 
for the small bias in Ai from charm decays on the -Stag 
side. The third distribution is centered at zero with a 
fixed width of 8ps. We have verified in simulation that 
the parameters of TZ{5t,aAt) for _B" — > K';lTr'^j events 
are similar to those obtained from the i?flav sample, even 
though the distributions of (TAt differ considerably. We 
therefore extract these parameters from a fit to the Bflav 
sample. We assume that the continuum background con- 
sists of prompt decays only and find that the At distri- 
bution is well described by a resolution function with the 
same functional form as used for signal events. The pa- 
rameters of the background function are determined in 
the fit. 

Figure shows the background-subtracted distribu- 
tions for mEs and AE for the selected B^ K*^^ can- 
didates. The background subtraction is performed with 
the event weighting technique described in |2^. Events 
contribute according to a weight constructed from the co- 
variance matrix for the signal, BB background and con- 
tinuum background yields and the probability Vs, V^g 
and Vqq for the event, computed without the use of the 
variable that is being displayed. The curves in the figure 
represent the signal PDFs used in the fit. FigureElshows 
the background-subtracted distributions of At for B"- 
and -B'^-tagged events, and the asymmetry as a function 
of At. 

In the fit to the B" K*'^^ sample we find 157 ± 16 
signal events, with 

Sk>o^ = -0.21 ±0.40 ±0.05 

and 

Ck-o^ = -0.40 ±0.23 ±0.03, 




nigg [GeV/c^] AE[GeV] 



FIG. 1: Background-subtracted distribution for rriES (left) 
and A.E (right) for 0.8 < mj^o^^o < 1.0. The lines represent 
the signal PDFs used in the fit, normalized to the fitted yield. 
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FIG. 2: Background-subtracted distribution for At with -Btag 
tagged as i?" (top) or (center), and the asymmetry 
^^0^0 (At) (bottom). The curves represent the PDFs for 
signal decays in the likelihood fit, normalized to the final fit 
result. 

where the first error is statistical and the second system- 
atic. The systematic uncertainties are described below. 
The extracted BB background is 9 ± 13 events. The lin- 
ear correlation coefficient between Sk'O^ and C/fo^ is 
0.07. The value of Cj^'O-y is consistent with the expec- 
tation of no direct CP violation. Since its uncertainty 
is much larger than that obtained from the partial rate 
asymmetry in self-tagging decays jl4j . we also perform 
the fit with Ck'Oj fixed to zero and find 

SK-o-yiC = 0) = -0.22 ± 0.42 ± 0.05. 

The counterintuitive increase in the error on Sx*Oj is a 
consequence of the likelihood contours in the S-C plane, 
shown in figure |31 not being perfectly ellipsoidal. 

Figure 0] shows the background-subtracted K'^tt'^ in- 
variant mass distribution for B^ K'^tt^j candidates. 
The K* (892) resonance is clearly visible and there is some 
evidence for the (1430). Figure |31 shows the back- 
ground subtracted distributions for tties and AE events 
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FIG. 5: Background-subtracted distribution for rriES (left) 
and Ai? (right) for 1.1 < m^o^o < 1.8. The hnes represent 
the signal PDFs used in the fit, normalized to the fitted yield. 



FIG. 3: Constant-likelihood contours in the S-C plane for 
^0 _^ ^*o^ corresponding to -2AlogL = 1, 2 and 3. The 
dashed circle is the physical boundary. 
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FIG. 4: Background-subtracted distribution for m^^o^o. For 
events with nij^o^o > 1.1 GeV/c^ the cut on the cosine of the 
helicity angle cos9k' is not applied. 



in the range 1.1 < toj<-o^o < 1.8 GeV/c^. In the fit to 
this sample we find 59 ± 13 signal events with 



S 



and 



= 0.9 ± 1.0 ±0.2 



-1.0±0.5±0.2, 



and 130 ±40 BB background events. The linear correla- 
tion coefficient between S^o-nO-y and Cxo-n^-y is —0.09. 

We consider several sources of systematic uncertainties 
related to the level and possible asymmetry of the back- 
ground contribution from other B decays. We evaluate 
this contribution using simulated samples of generic B 
decays and of generic B Xs^ decays. For the latter 
we use the Kagan-Neubert model for the photon en- 
ergy spectrum and JETSET for the fragmentation of the 
s quark. Since the final state multiplicity predicted by 
the fragmentation model is significantly different from a 



recent BABAR measurement [2j| , we reweight events ac- 
cording to their multiplicity. Prom these studies we esti- 
mate about 30 (140) events in the K* (non-iiT*) sample, 
with approximately equal contributions B — > Xg'y decays 
and other (generic) B decays. In the K* sample we fit a 
contribution from BB background consistent both with 
zero and with the rate predicted from the simulation. In 
the non-i^T* sample the fitted contribution from other B 
decays is significantly larger, as expected from the sim- 
ulation. Using the Monte Carlo estimates for the yields, 
we assess the impact of a potential CPV asymmetry in 
the BB background by varying Sg^ and C^g within an 
appropriate range that is derived from the composition of 
the BB background sample. We assign a systematic un- 
certainty of 0.04 (0.03) on S (C) in the K* sample and 
an uncertainty of 0.2 for both S and C in the non-ii'* 
sample. 

We quantify possible systematic effects due to the ver- 
tex reconstruction method in the same manner as in 
Ref. [T9|. estimating systematic uncertainties on S (C) 
of 0.023 (0.014) due to the vertex reconstruction tech- 
nique and uncertainties in the resolution function, and 
0.020 (0.007) due to possible misahgnments of the SVT. 
Finally, we include a systematic uncertainty due to im- 
perfect knowledge of the PDFs used in the fit, which 
amounts to 0.02 (0.01) for the K* (non-K*) sample. 

In summary, we have performed a new measurement 
of the time-dependent CPV asymmetry in S° K*'^j 
decays. Within the large statistical uncertainties our 
measurement is consistent with the SM expectation of a 
small CPV asymmetry and with other measurements . 
We have also explored the possibility of measuring the 
CPV asymmetry in the region with a K'^iir'^ invariant 
mass above the K*'^ region, 1.1 < nii^a^o < 1.8 GeV/c^. 
We find that the signal yield, though consistent with 
the expectation, is too small for a meaningful asymme- 
try measurement. These results supersede our previous 
measurement which was based on a subset of the data 
presented here. 
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